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Abstract 

Observations show that small-amplitude prominence oscillations are usually damped 
after a few periods. This phenomenon has been theoretically investigated in terms 
of non-ideal magnetoacoustic waves, non-adiabatic effects being the best candidates 
to explain the damping in the case of slow modes. We study the attenuation of 
non-adiabatic magnetoacoustic waves in a slab prominence embedded in the coro- 
nal medium. We assume an equilibrium configuration with a transverse magnetic 
field to the slab axis and investigate wave damping by thermal conduction and 
radiative losses. The magnetohydrodynamic equations are considered in their lin- 
earised form and terms representing thermal conduction, radiation and heating are 
included in the energy equation. The differential equations that govern linear slow 
and fast modes are numerically solved to obtain the complex oscillatory frequency 
and the corresponding eigenfunctions. We find that coronal thermal conduction and 
radiative losses from the prominence plasma reveal as the most relevant damping 
mechanisms. Both mechanisms govern together the attenuation of hybrid modes, 
whereas prominence radiation is responsible for the damping of internal modes and 
coronal conduction essentially dominates the attenuation of external modes. In addi- 
tion, the energy transfer between the prominence and the corona caused by thermal 
conduction has a noticeable effect on the wave stability, radiative losses from the 
prominence plasma being of paramount importance for the thermal stability of fast 
modes. We conclude that slow modes are efficiently damped, with damping times 
compatible with observations. On the contrary, fast modes are less attenuated by 
non-adiabatic effects and their damping times are several orders of magnitude larger 
than those observed. The presence of the corona causes a decrease of the damping 
times with respect to those of an isolated prominence slab, but its effect is still 
insufficient to obtain damping times of the order of the period in the case of fast 
modes. 
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i 1 Introduction 
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Solar prominences are large-scale coronal magnetic structures whose material, 
cooler and denser than the typical coronal medium, is in plasma state. Promi- 
nences are supported against gravity by the coronal magnetic field, which also 
maintains the prominence material thermally isolated from the corona. Small- 
amplitude osci llations in solar prominences were detected almost 40 years ago 
( iHarveyi . Il969l ). These oscillatory motions seem to be of local nature and their 
velocity amplitude is typically less than 2-3 km s _1 . Observations have also al- 



lowed to measure a wide range of periods between 30 s (IBalthasar et all Il993f ) 



and 12 h (IFoullon et all 12004 ) . More recently, some high-resolution observa- 



tions of p rominence oscil l ations by the Hino d e/SOT instrument have bee n 



reported flOkamoto et all . 120071 : iBerger et all 120081 : lOfman fc Wand . 120081 ) 



From the theoretical point of view, the oscillations have been interpreted by 
means of the magnetoacoustic ei genmodes suppo r ted by the prominence body. 
A rece nt example is the work by iTerradas et al.l (120081 ) in which the observa- 
tions of lOkamoto et al.l (120071) are int e rpreted as fast k i nk waves. The reader is 
referred to Toiiver k Ballesterl (l2002f ): iBallesterl (l2006h : iBaneriee et~aD (120071 ) 
for extensive reviews of both observational and theoretical studies. 
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Evidence of the attenuation of small-amplitude prom i nence oscillations has 



been reported in some works (IMolowny-Horas et al.l . Il999l ; ITerradas et al. 



20021 ; iLinl . 120041 ). A typical feature of these observations is that the oscillatory 
motions disappear after a few periods, hence they are quickly damped by one 
or several mechanisms. The theoretical investigation of this phenomenon in 
terms of magnetohydrodynamic (MHD) waves has been broached by some au- 
thors by removing the ideal assumption and by including dissipative terms in 
the basic equations. Non-adiabatic effects appear to be very efficient damping 
mechani sms and have been investiga t ed with the help of simple prom inence 
models dBallail . 120031 : ICarbonell etall 12004! . 120061 : ITerradas et aD . l2005h . Nev- 
erth eless, other dampin g mechanisms have been al so proposed, like wave leak- 
age (ISchuteensl. Il997alib ; ISchutgens fc Tothl . Il999l ). dissi pation by ion-neutra l 
collisions (jForteza et al. . 2007 ) and resonant absorption ( Arregui et al|[2008l ). 



In a previous work (jSoler et al.l . 120071 . hereafter Paper I), we have studied for 
the first time the wave attenuation by non-adiabatic effects of a prominence 
slab embedded in the corona. In that work the magnetic field is parallel to 
the slab axis and it is found that the corona has no influence on the internal 
slow modes, but it is of paramount importance to explain the damping of fast 
modes, which are more attenuated than in simple models that do not consider 
the coronal medium. Following the path initiated in Paper I, here we investi- 
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gate the wave damping due to non-adiabatic mechanisms (radiative losses and 
thermal conduction) in an equilibrium made of a prominence slab embedded 
in a coronal medium, but now we consider a magnetic field transverse to the 
slab axis. This configuration and that studied in Paper I correspond to limit 
cases, since measurements with Zeeman and Hanle effects indicate that the 
magnetic field lines are skewed to the long axis of prominences. On average, 
the prominence axis and the magnetic field form an angle of about 20 deg. 
Thus, the skewed case is relegated to a future investigation. 



The eq uilib rium configuration a ssumed here was analysed in detail by lJoarder fc Roberts 
(119921 ) and lOliver et al.l (119931 ) in the case of ideal, adiabatic pertu rbations. 



The m ain difference between both works is in the treatment of gravity. I Joarder fc Roberts 
(119921 ) neglected the effec t of gravity and so straight field lines were consid- 
ered. On the other hand, lOliver et al.l (119931 ) took gravity into account and 
assumed curved field lines according to the Kippenhahn &: Schiilter ( 1957t) 
mode l modified to include the surrounding coronal plasma (IPoland fc Anzerl . 
19711 ) . Despite this difference, both studies agree in establishing a distinction 
between different normal modes depending on the dominant medium sup- 
porting the oscillation. Hence, internal modes are essentially supported by 
the prominence slab whereas external modes arise from the presence of the 
corona. In addition, hybrid (or string) modes appear due to the combined 
effect of both media. 
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The investigation of the thermal attenuation of oscillations supported by such 
equilibrium is unsettled to date and, in deed, this is the main motiv ation for 



the present study. However, two works (ITerradas et all l200ll. 120051) st udied 



the wave damping in an isolated prominence slab. ITerradas et al.l (1200 ll ) con- 
sidered radiative losses given by the Newton ian law of cooling as damping 



mec hanism a nd stu died the attenuation in the iKippenhahn fc Schiilter 



and iMenzell (ll95ll ) prominence models. Subsequently, ITerradas et al 



considered a more complete energy equation including optically thin radiation, 
plasma heating and parallel thermal conduction, and assumed straight field 
lines since gravity was neglected. The main conclusion of both works is that 
non-adiabatic mechanisms are only efficient in damping slow modes whereas 
fast modes remain almost undamped. Nevertheless, in the light of the results 
of Paper I, the presence of the coronal medium can have an important reper- 
cussion on the wave damping. The investigation of this e ffect is the main aim o f 
the present work. Therefore, we extend here the work of ITerradas et al.l (120051 ) 
by considering the presence of the corona and neglect the effect of gravity as 
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Joarder fc Roberts! (119921 ) for simplicity. 



This paper is organised as follows. Section [2] contains a description of the 
equilibrium configuration and the basic equations which govern non-adiabatic 
magnetoacoustic waves. Then, the results of this work are extensively discussed 
in Sect. EJ Finally, our conclusions are given in Sect. |U 
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8i 2 Equilibrium and basic equations 



The equilibrium configuration (see Fig. [TJ is made of a homogeneous plasma 
slab with prominence conditions (density p v and temperature T p ), whose axis is 
orientated along the z-direction, embedded in a coronal environment (density 
p c and temperature T c ). The system is bounded in the x-direction due to 
the presence of two rigid walls representing the solar photosphere, but it is 
unlimited in the y- and z-directions. The width of the prominence slab is 2x p 
and the total width of the system is 2x c . The magnetic field is transverse to 
the prominence slab, B = B e x , with B everywhere constant. 
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Fig. 1. Sketch of the equilibrium. The dark region represents the prominence slab 
while the light region corresponds to the corona. The photospheric walls are the two 
hatched areas on both sides of the corona. 

90 In order to find the basic equations tha t govern non-ad i abatic magnetoacoustic 

91 waves we follow the same process as in lTerradas et al.l (120051) . We consider the 

92 usual MHD equations (Eqs. (l)-(6) of iTerradas et all 120051 ) in which non- 
93 adiabatic terms have been included in the energy equation, 



Dp _ jpDp. 

Dt p Dt 17 



l)[pL(p,T)-V-(K-VT)} = 0, 



where p, p and T are the gas pressure, density and temperature, respectively. 
The quantity 7 is the adiabatic ratio, her e taken 7 = 5/ 3 . The non-ideal 
terms in Eq. ([TJ are explained in detail in ICarbonell et al.l (120041 ) . Thermal 
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conduction is represented by V • (k • VT), where k is the conductivity tensor 
which in coronal and prominence applications is usually approximated by its 
parallel component to the mag netic field, k\\ = 10~ 11 T 5 / 2 W ur 1 K -1 . Radia- 
tive losses and heating are evaluated together through the heat-loss function, 
L(p, T) = x*pT a — hp a T b , where radiation is parametrised with x* an d ot (see 
Table I of Paper I) and the heating scenario is given by exponents a and h 
flRosner et all Il978l : iDahlburg fc Mariskal . Il988h . 



105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 



Regarding our equilibrium configuration, the reader must be aware that, al- 
though there have been some attempts to construct a self-consisten t promi- 



nence model including both magnetostatics and thermodynamics (e.g. lMilne at al. 



19791 ; iLow fc Wul . Il98ll ; lAnzer fc Heinzell . Il999l ). to date this task remains to 



be done. Here, we consider a simplified prominence-corona configuration, but 
it includes the two basic ingredients observed in real prominences. First, the 
existence of a steep temperature gradient between the prominence and the 
corona and, second, the apparent thermal isolation of the prominence mate- 
rial from the much hotter corona. The first point is addressed by considering 
that the temperature profile is a step function, and so the prominence- corona 
transition region (PCTR) has n ot been considered. This choice is supported 
by results of previous works (e.g. lOliver fc Ballesterl . Il996l ) which showed that 
the PCTR has a minor influence on the prominence oscillatory modes. On 
the other hand, to represent the thermal isolation we have neglected the heat 
flux due to thermal conduction at the boundary between the prominence and 
the corona. Therefore, we impose that both the prominence and the corona 
are isothermal and thermally isolated, and so radiative losses and heating are 
locally balanced, i.e. L(p ,T ) = 0, where p and T are the local equilibrium 
density and temperature, respectively. 



124 Assuming that the plasma is at rest in the equilibrium state (i.e. no flux of 

125 material) and considering small pertur bations, we find the li nearised version of 

126 the MHD equations (Eqs. (10)-(15) of lTerradas et all 120051 ). According to the 

127 geometry of our model, we assume perturbations of the form f±(x) expi(ut + 

128 k y y + k z z) and exclude Alfven waves from this analysis by considering only 

129 motions and propagation in the xz-plane (v y = 0, k y = 0). Now we combine 

130 the resultant expressions and eliminate all the perturbed quantities in favour 

131 of the velocity perturbations, v x and v z , and the temperature perturbation, 

132 T\. By this process, we obtain three coupled ordinary differential equations, 



d 2 ?; 

c s 77 + 1^ v x + ik z c s 
dx 2 dx 



2 dv z 



iuic^ dTi 
T dx 



(2) 



V A 



2 d 2 v z 



dx 2 



+ 



2 ; 2 




+ i k C * dV * 



, c 2 s T x 
, +uk z —— 
7 dx 7 T 



0. 



(3) 
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0, (4) 



2 



135 where = ^ is the adiabatic sound speed squared whereas v\ = — is 

136 the Alfven speed squared, po and B denote the equilibrium gas pressure 

137 and magnetic field strength, respectively, and /i is the magnetic permittivity 

138 (n = 47rl0 -7 in MKS units). Quantities ojt and u p are defined as follows, 

139 Up = — (L + PqL p ) , w T = —T L T , 

Po Po 



140 L p , Lt being the partial derivatives of the heat-loss function with respect to 

141 density and temperature, respectively, 



Lt 
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Equations (j2J), (0]) and (TjJ govern fast and slow magnetoacoustic waves to- 
gether with the thermal or condensation mode. In this work we do not study 
the thermal wave sin ce we pay our attention to the magnetoacoustic modes. 



Terradas et al.l (120051 ) found an approximate analytical solution of Eqs. 



by neglecting thermal conduction, a valid assumption in prominence plasmas. 
However, thermal conduction has an important role in coronal conditions and 
cannot be neglected in order to perform a realistic description of the oscil- 
latory modes supported by our equilibrium configuration. He nce, w e solve 
the full set of Eqs. ©-Q u sing the numer i cal co de PDE2D flSewelll . l2003h 



based on finite elements (see iTerradas et al.l . 120051 . for an explanation of the 
method). The jump conditions at the interface between the prominence and 
the corona are automatically we ll-treated by the code. These jump conditions 
are flGoedbloed fc Poedtsl . Eopj ) : 



0. 



B 



0, [p] = 0, 



(5) 



157 where v and B are the perturbed velocity and the magnetic field vectors, re- 

158 spectively. For a complete closure of the system we need to supply a physically 

159 consistent set of boundary conditions for the perturbations at the photospheric 
wo walls, x = ±x c . In this work, we consider two different sets of boundary con- 

161 ditions, 

162 Vx = Vz = Tx = 0, at x = ±x c , (6) 
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and 



T[ = 0, at x = ±x c 



(7) 



165 where ' indicates derivative with respect to x. Both sets consider line-tied 

lee conditions for the velocity perturbations, i.e. the disturbances are unable to 

167 perturb the dense photospheric plasma which acts as perfectly rigid wall. On 

168 the other hand, sets © and (J7J) differ by the condition for T\, which has 

169 different physical implications. Set ([6]) assumes that the perturbation to the 

170 temperature vanishes at x — ±x c and this means that the photospheric walls 

171 are taken as isothermal. On the contrary, set ([7j) considers a zero-temperature 

172 gradient for the perturbation between the corona and the photosphere, so no 

173 perturbed heat flux is allowed at the boundaries. From our point of view, 

174 set <Q makes more physical sense than set ([7]), since one can expect that the 

175 much denser photospheric plasma can instantaneously radiate away any in- 

176 coming perturbed heat flux from the corona. However, set ([7j) imposes that 

177 there is no heat exchange between the corona and the photosphere, although 
the temperature perturbation c an have a non-ze r o valu e at the walls. Regard- 



178 



180 



179 ing these boundary conditions, ICargill fc Hoodl (119891 ) performed a study of 



the thermal stability of wave and thermal modes in a Cartesian coronal slab 

181 and pointed out that the solutions computed by assuming the boundary con- 

182 ditions given by set ([H]) are more thermally stable than those obtained for 

183 boundary conditions of set ([7]). 

184 For a fixed real k z , the numerical solution of Eqs. (T5]), ([3]) and (@J provides 

185 with a complex frequency, to = uj^ + iui. In the ideal, adiab atic case ui-\ = and 



186 t heref ore the solutions of Eqs. (T5J), ([3]) and (jl]) are those of iJoarder fc Roberts 



187 (119921 ). Using the real and imaginary parts of the frequency, we can com- 

188 pute the oscillatory period, P, the damping time, td, and the ratio of both 

189 quantities, 

uo P = — , r D = — , — = . (8) 

Co>r UJ\ P Z7T 0J\ 



i9i 3 Results 



192 Unless otherwise stated, the following equilibrium parameters are considered in 

193 all computations: T p = 8000 K, p p = 5 x 10~ n kg m" 3 , T c = 10 6 K, p c = 2.5 x 

194 10 -13 kg m~ 3 , Bq = 5 G, x p = 3000 km and x c = 10x p . The coronal density is 

195 computed by fixing the coronal temperature and imposing pressure continuity 
we across the interfaces. In addition, we assume an optically thin prominence 

197 plasma (regime Prominence (1) of Paper I) and a constant heating per unit 

198 volume (a = b = 0). In all the following expressions, subscript indicates 
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199 local equilibrium values, while subscripts p and c denote quantities explicitly 

200 computed with prominence and coronal parameters, respectively. 



201 3.1 Dispersion diagram and wave modes 



202 
203 
204 
205 
206 
207 
208 
209 
210 
211 



Solutions of Eqs. (j2j)-(jlj) can be grouped in internal, external and hybrid 
modes. Although there is an infinite number of harmonics for internal and 
external modes, only two hybrid modes are possible: th e hybrid slow mod e 
and the hybrid fast mode (this nomenclature is taken from lOliver et all 119931 ). 



Figure [2] shows the dimensionless real part of the frequency versus k z x p for 
the fundamental symmetric oscillatory modes (i.e. solutions with v x even with 
respect x — 0) and some of their harmonics, where we have assumed the 
boundary conditions given by Eq. ([6]). A similar diagram can be obtained for 
the antisymmetric modes (i.e. solutions with v x odd with respect x — 0) and 
for the other set of boundary conditions (Eq. (j7|)). 



30 



20 



10 



10 



Fig. 2. Dimensionless real part of the frequency versus k z x p for the oscillatory 
symmetric modes: hybrid slow (solid line at the bottom), fundamental internal slow 
and first harmonics (dotted lines), fundamental external slow and first harmonics 
(dashed lines), fundamental internal fast and first harmonic (dash-dotted lines) and 
fundamental external fast (three dot-dashed line at the upper left corner). 

Th e behaviour of the real part of the frequency is the same as that explained 
by lOliver et al.l (I1993I ) . The value of for both internal and external slow 
modes and for the hybrid slow mode shows a very weak dependence on k z 
since almost horizontal lines are seen in Fig. [2J On the contrary, both internal 
and external fast modes show a quasi-parabolic dependence on k z (this is 
also applicable to the hybrid fast mode, present in the dispersion di agram for 
antisymmetric modes). The reader is referred to lOliver et al.l (I1993I ) for more 
extensive details about the behaviour of cur. 
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Hybrid slow mode 




Fig. 3. Modulus of the eigenfunctions v x , v z and T\ (in arbitrary units) versus the 
dimensionless distance to the slab axis corresponding to the fundamental oscilla- 
tory modes for k z x p = 1. The solid line corresponds to the boundary conditions 
v x = v z = T[ = at x = ±x c , while the dotted line corresponds to the boundary 
conditions v x = v z = T\ = at x = ±x c . The shaded region shows the location of 
the prominence slab. 



220 Next, we focus on the fundamental modes and their eigenfunctions v x , v z and 

221 Ti are displayed in Fig. [3] for k z x p = 1 and for both sets of boundary con- 

222 dit ions. The spatial st ructure of the disturbances v x and v z is the one shown 

223 by lOliver et al.l (119931 ). Hence, non-adiabatic effects do not modify the spa- 

224 tial behaviour of velocity perturbations. Internal modes produce large plasma 
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225 displacements inside the slab, external modes achieve large amplitudes in the 

226 corona and the amplitude of hybrid modes is of the same order in both media. 

227 It is worth to mention that the hybrid fast mode can be considered as an 

228 internal-like mode for large k z x p since the amplitude of its perturbations in 

229 the corona decreases as k z x p increases. Regarding the temperature perturba- 

230 tion, it is larger for the slow modes than for the fast modes and, in general, is 

231 larger in the prominence than in the corona. Finally, the differences between 

232 the eigenf unctions for the two sets of boundary conditions (Eqs.(JB])-(j7])) are 

233 only relevant for the hybrid slow mode. 

234 From the observational point of view, internal and hybrid modes could be more 

235 easily observed than external modes by instruments focusing on prominences, 

236 since the amplitude of the latter ones is very small in the prominence body. 

237 For this reason, the results corresponding to internal and hybrid modes are 

238 the most interesting for prominence seismology. However, here we study the 

239 three kinds of solutions in order to perform a complete description of the 

240 fundamental wave modes supported by the equilibrium configuration. 



24i 3.2 Mode coupling 



242 



Oliver et~aD dl993h showed that avoided crossings occur in the dispersion dia- 

243 gram when two solutions couple and interchange their magnetoacoustic prop- 

244 erties. Nevertheless, no avoided crossings seem to take place in our dispersion 

245 diagram (Fig. [2]) since the curves of cur for the internal fast modes and for 

246 the slow modes cut each other. This fact can be understood by considering 

247 that in the present, non-adiabatic case the complete dispersion diagram is in 

248 a three-dimensional space because the frequency has an imaginary part. So, 

249 Fig. factually corresponds to a projection of the complete three-dimensional 

250 dispersion diagram on the fc^R-plane. 

251 Upon exploring the complete dispersion diagram, we have found that three 

252 different couplings can take place: 

253 (1) If the imaginary parts of the frequency of the coupling modes differ by 

254 several orders of magnitude, there is no avoided crossing between the real 

255 parts. Hence, the coupling between modes is "weak" and only becomes 

256 apparent by means of a slight mutual approach of the imaginary parts of 

257 uj (see Fig. |U left panel). 

258 (2) If both imaginary parts of the frequency have a similar value, the real 

259 parts show an avoided crossing and so a "strong" coupling takes place 

260 (see Fig. HI mid panel). 

261 (3) In very peculiar cases, an "anomalous" coupling takes place when the 

262 imaginary parts of u of the two coupling modes repel each other (see 
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263 Fig. HI right panel). This situation has important effects on the wave 

264 stability, as we explain in Sect. 13.51 



265 The b ehaviour of the mode coupling was previously described by lTerradas et al 



266 ( 120011 ) in the cases th at we call "weak" and "strong" couplings (compare our 



Fig. 1 with Fig. 12 of lTerradas et all l200lh 






Fig. 4. Three-dimensional dispersion diagrams (solid lines) close to a coupling be- 
tween a fast mode and a slow mode. Dashed and dotted lines are the projections 
of the dispersion curves on the horizontal and vertical planes. The left-hand side 
panel presents a "weak" coupling, the middle panel shows a "strong" coupling and 
the right-hand side panel displays an "anomalous" coupling. 



268 3.3 Periods and damping times 
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Hereafter we restrict ourselves to the fundamental modes and compute the 
oscillatory period, P, the damping time, r D , and the ratio of the damping 
time to the period as functions of the dimensionless wavenumber, k z x p . We 
consider values for k z x p between 0.01 and 3, which correspond to wavelengths 
between 5 x 10 3 km and 10 5 km, approximately. These valu es cover the range of 
typic ally observed wavelengths in prominence oscillations (lOliver &: Ballesterl . 
20021 ) . The results of the computations are displayed in Fig. [5] considering the 
two sets of boundary conditions (Eqs. (EJ)-(J71)). 

The periods obtained he re a gree with those pro vided and commented by 
Joarder fc Roberts! (119921 ) and lOliver et al.l (119931 ) , therefore we turn our at- 
tention to the damping times. Regarding slow modes, we see that they are 
strongly damped, with values of td/P close to 1 for the three modes. How- 
ever, fast modes are much less attenuated and the obtained values of t^/P are 
much larger than those observed. This fact involves an important difference 
with the results of Paper I, in which fast waves were efficiently attenuated 
for some values of the wavenumber. Such as happens with the period, the 
damping time of slow modes is almost independent of k z . On the contrary, 
fast modes are less attenuated for small k z than for large k z . This evidence 
can be understood by means of the following arguments. Considering k z = 



11 



£ no 



105 - 
100 _ 
0.0 



15 
1? 10 



Hybrid slow mode 

12001 1 1 ' ' ' 



0.5 1.0 1.5 2.0 2.5 3.0 

k, X, 



s 1000 
j 800 
- 600 

400 

200 E 

0.0 0.5 /.0 r.5 2.0 2.5 3.0 

Hybrid fast mode 

10000? 



,g 1000 



10 f ■ ■ ~ ■ ■ ; 

s • 

G - 

4 - z 

2 - : 

[ , , , , , : 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 
k, x. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



25 r~ 
24 - 
23 K. 
22 - 
21 - 
20 _ 
0.0 



0.5 1.0 1.5 2.0 2.5 3.0 





4 




3 








2 


a. 


1 








0.5 1.0 1.5 2.0 2.5 3.0 



-s 8 
C 



0.5 1.0 1.5 2.0 2.5 3.0 
k. x. 





Fund, internal slow mode 

75 



_ 70 
| 05 
^ 60 

56 

50 

0.0 0.5 >.0 1.5 2.0 2.5 3.0 

Fund, internal fast mode 

10 s 




Fund, external slow mode 

12 



0.0 0.5 1.0 1.5 2.0 2.5 3.1 
k z x 

Fund, external fast mode 
10' 

? 

J. io l 

° 10' 
10' 



0.0 0.5 1.0 1.5 2.0 2.5 3.1 

K x. 



0.5 1.0 1.5 2.0 2.5 3 

k, x. 



^ 3 ' 
t? 2 - 



0.0 0.5 1.0 1.5 2.0 2.5 .3.0 
k. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
k, i„ 



0° 
0" 
s 
0" 
0' 

0.0 0.5 1.0 1.5 2 2.5 3.0 
k x 



Fig. 5. Period (left), damping time (centre) and ratio of the damping time to the 
period (right) versus k z x p for the fundamental oscillatory modes. The solid line 
corresponds to the boundary conditions v x = v z = T[ = 0, while the dotted line 
corresponds to the boundary conditions v x = v z = T\ = 0. 

2B8 then Eqs. (J2J) — become, 

9 d 2 t> x 9 iujcl dTi , , 

" C ^ + ^-^d^^°' (9) 

vl^+u 2 v z = 0, (10) 
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292 Equations © and (TTTT) are still coupled and govern slow and thermal waves, 

293 which are affected by non-adiabatic mechanisms through the terms with k», 

294 UJj' and LJ p in Eq. fill I) . On the contrary, Eq. (jTUJ) is now decoupled from the 

295 rest and governs fast modes alone, which become pure Alfven waves and are 

296 not affected by non-adiabatic terms. Thus, for k z — > fast waves tend to 

297 the ideal, undamped behaviour. When k z is increased, fast modes are more 

298 affected by acoustic effects and their damping time decreases and stabilises. 

299 The little peaks shown in the bottom panels of Fig. corresponding to the 

300 external fast mode, are in fact the result of "strong" couplings with slow mode 

301 harmonics. The differences arising from the different boundary conditions are 

302 only of importance for the hybrid slow mode, as we indicated in Sect. 13.11 

303 We see that the boundary condition T[ = produces a substantially stronger 

304 damping for the hybrid slow mode than the condition T\ = 0. 

305 Finally, an approximate value to the frequency of internal and external slow 

306 modes can be obtained by considering the approximation given in App. B of 

307 Paper I, namely 

308 to ~ Ak x , (12) 



309 where k x is the wavenumber in the field direction and A is the modified sound 

310 speed due to the presence of non-adiabatic effects, defined in Paper I as follows 

■( 7 -i) (a«n** + 



A = 



7 



(7-1) {^kl + UT)+iu 



(13) 



The value of k x is fixed by the equilibrium geometry, but for simplicity we 
consider now the analytical approx imations of the dominant wavenumbers 
given by IJoarder fc Roberts! (119921 ) in the adiabatic case and for the long 
wavelength limit, namely 



h, i 



71 



(14) 



317 for the fundamental internal mode, and 

71 

318 k x ps (15) 



319 for the fundamental external mode. One must bear in mind that the wavenum- 

320 bers in the present, non-adiabatic case are complex quantities, but we expect 
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321 that their real part is similar to that in the adiabatic case, such as happens 

322 with the value of the frequency. Applying now Eq. (Tl2|) to the internal slow 

323 mode, i.e. considering prominence parameters in the expression for A (Eq. IT5T) 

324 and the approximation for k x given by Eq. (1141) . one obtains P ~ 23.50 min, 

325 td ~ 73.05 min and td/P ~ 3.11. On the other hand, if the process is re- 

326 peated for the external slow mode, this gives P ~ 6.20 min, td ~ 6.70 min 

327 and tb/P ~ 1.08. We see that these approximate values reasonably agree with 

328 those numerically obtained and represented in Fig. [51 



329 3.4 Importance of the damping mechanisms 




Fig. 6. Damping time versus k z x p for the fundamental modes: a) hybrid slow, b) 
internal slow, c) external slow, d) hybrid fast, e) internal fast, and /) external fast. 
Different linestyles represent the omitted mechanism: all mechanisms considered 
(solid line), prominence conduction eliminated (dotted line), prominence radiation 
eliminated (dashed line), coronal conduction eliminated (dot-dashed line) and coro- 
nal radiation eliminated (three dot-dashed line). Arrows in panels d and e point the 
location of thermal instabilities {uj\ < 0) which appear if prominence radiation is 
omitted (dashed line). 

330 In order to know which are the mechanisms responsible for the damping of 

331 each mode, we now follow the same procedure as in Paper I. We compare the 

332 damping time obtained when considering all non-adiabatic terms (displayed in 

333 the middle column of Fig. [5]) with the results obtained when a specific mecha- 

334 nism is removed from the energy equation (Eq. (CQ)). This analysis allows us to 

335 know whether the omitted mechanism has a relevant effect on the attenuation. 

336 Before undertaking this investigation, we need to know if both sets of boundary 

337 conditions are adequate in the absence of thermal conduction. If one imposes 
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338 km = in Eq. (j4j) then T\ can be written as function of v z and v' x 



x • 



T {l + iu p /uj) , 
Tl = , • ft — TT + lkzVz > ' ( 16 ) 



340 which can be substituted into Eqs. (j2j) and fl3j) in order to obtain two coupled 

341 differential equations involving the perturbed velocities alone, 

9 / iuj\ d 2 v x 9 , 9 / ^\ df 2 . . 

- c H 1 + ^)d^ + ^ + * c -( 1 + ^)d7 = °- (17) 



dV 
dx 2 



2 Cg / IUO 
«A+-(l + — 

7 V 

ZOA dl^ 



+ 1+ , , , 

7 V / ax 



0. 



(18) 



343 Here Cj is introduced to simplify the notation, 
uj t + ia»/ (7 - 1) 



UJ 



1 + iuj p /uj 



345 Now, the system formed by Eqs. (TTTT) and (TTHT) is fully determined by assuming 

346 only boundary conditions for v x and v z . Hence, the behaviour of 7\ at the 

347 boundaries cannot be imposed but is fixed by the conditions over the velocity 

348 perturbations. If one takes v x = v z = T[ = as boundary conditions, then 

349 Eq. (1TB]) yields the constraint v" x + ik z v' z = 0, which substituted in Eq. ( TTTT) 

350 automatically gives the redundant condition v x = 0. On the other hand, if 

351 one assumes v x = v z = T\ = at x = ±x c , then Eq. (TT6|) now imposes 

352 v' x = at the boundaries. This last condition substituted in Eq. (TTHT) gives 

353 the extra condition v" z = over the system, which implies a new restriction 

354 that is not generally satisfied by all solutions. Thus, T[ = reveals itself 

355 as the "natural" boundary condition for the temperature perturbation when 

356 thermal conduction is neglected. So, for the following investigation we restrict 

357 ourselves to the boundary conditions v x = v z = T[ = since the conditions 

358 v x = v z = T\ = are not consistent with the differential equations when 

359 thermal conduction is neglected. 

360 The results of the computations are displayed in Fig. El Although we have 

361 explored a wide range of values of k z , the plots are only drawn again for 

362 0.01 < k z x p < 3 since we have found that the importance of the damping 

363 mechanisms does not show a strong dependence on k z . Regarding slow modes, 

364 we clearly see that the damping of the internal mode is dominated by the radi- 

365 ation from the prominence plasma, as expected, while coronal conduction has 
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366 a minor effect. On the other hand, the hybrid and external modes are affected 

367 by coronal conduction together with prominence radiation. Both mechanisms 

368 have a similar influence on the hybrid mode, while coronal conduction domi- 

369 nates the attenuation of the external mode. This result for the hybrid mode 

370 is coherent with the fact that its perturbations achieve large amplitudes both 

371 in the prominence and the corona (see top row of Fig. [3]), so one expects that 

372 the most relevant damping mechanisms of each medium govern together the 

373 attenuation of the hybrid mode. However, the result for the external mode is 

374 a priori surprising because its perturbations are very small in the prominence 

375 (see fifth row of Fig. [3]) and one expects that the prominence-related mecha- 

376 nisms have a minor effect on its damping. The following discussion attempts 

377 to explain why prominence radiation affects so much the external mode. 

378 The equilibrium configuration assumed in the present work implies an addi- 

379 tional complication with respect to the equilibrium considered in Paper I, in 

380 which magnetic field lines were taken parallel to the interface between the 

381 prominence and the corona. Hence, both media were thermally isolated in the 

382 model of Paper I since there was no transfer of energy from one medium to the 

383 other. However, in the present model thermal conduction connects both media 

384 since field lines are transverse to the interfaces. This fact allows heat transfer 

385 between the prominence and the corona. So, some energy can flow along field 

386 lines and can be injected from the corona into the prominence, where the en- 

387 ergy is efficiently radiated away by the plasma. In this way, the influence of 

388 prominence radiation on the damping of the external slow mode, and also the 

389 hybrid slow mode, is amplified by means of coronal thermal conduction. 

390 Next we turn our attention to the fast modes. At first sight, the behaviour of 

391 the fast modes when a specific mechanism is removed from the energy equation 

392 is absolutely different from that seen in the case of the slow modes and needs 

393 more extensive explanations. In Sect. 13.3} we commented that the damping 

394 time of the fast modes is affected by the couplings with the slow modes. 

395 Now, we see that the nature of these couplings (being "weak", "strong" or 

396 "anomalous") changes depending on which is the non-adiabatic mechanism 

397 omitted in the energy equation. These changes in the coupling nature cause 

398 the damping time of the hybrid fast mode and the internal fast mode to vary 

399 from small values to very large values depending on the proximity to the 

400 couplings. So, we see that the consideration of both prominence radiation and 

401 coronal conduction has the effect of smoothing the curves of tq. 

402 In addition, the results corresponding to hybrid and internal fast modes show 

403 the appearance of thermal instabilities in very localised values of k z x p when 

404 prominence radiation is neglected (dashed lines), since then the interactions 

405 between fast modes and external slow modes leads to "anomalous" couplings. 

406 At these couplings, the value of uj\ for the fast modes is pushed towards neg- 

407 ative values (see the right-hand panel of Fig. 0J. Such a situation has very 
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408 important repercussions on the wave behaviour since for uo\ < waves are 

409 amplified in time. The location of these instabilities in panels d and e of Fig. [6] 

410 have been pointed by means of arrows. 



4ii 3.5 Wave instabilities 



412 Wave instabi lities discuss ed in Sect. 13.41 require a more in-depth investigation. 

413 According to lFieldl (119651 ). the criterion for the appearance of wave instabilities 

414 is given by 



K\\ 



Po 



-kl + L T + 



7 



(19) 



where k x is the wavenumber in the field direction. Results of ICarbonell et al. 
( 120041 ) . see also Paper I, point out that the heating scenario used in our cal- 



culations (constant heating per unit volume) cannot lead to thermal desta- 
bilisation. So, we can affirm that instabilities described in Sect. 13.41 are not 
caused by the heating mechanism. In addition, instabilities only appear when 
radiative losses are omitted. In such situation, the instability criterion becomes 



K 



—k 2 x < 0. 
Po 



(20) 



423 Equation (120]) is never satisfied unless an additional source of heating is 

424 present, which seems to be the present case. This extra energy source cor- 

425 responds to heat injected from the corona into the prominence by thermal 

426 conduction, as was commented in Sect. 13.41 In the absence of radiation, promi- 

427 nence thermal conduction is the only mechanism that can dissipate this extra 

428 injected heat. One expects that in such situation the value of k x grows in 

429 order to increase the efficiency of prominence conduction. Figure [7| shows the 

430 eigenfunction of the temperature perturbation corresponding to the internal 

431 fast mode for k z x v ~ 0.3 when all non-adiabatic mechanisms are considered, 

432 panel a) , and when radiative losses from the prominence plasma are omitted, 

433 panel b). For this value of k z x p , the wave becomes unstable (uj < 0) if promi- 

434 nence radiation is omitted. We see that smaller spatial-scales (i.e. larger k x ) 

435 are obtained within the prominence when prominence radiation is not taken 

436 into account, as expected. Although the efficiency of prominence conduction 

437 is increased in this way, it is still not enough to stabilise the perturbation. 

438 This last discussion points out that prominence radiative losses are of paramount 

439 importance to stabilise the disturbances. The efficiency of prominence radia- 

440 tion can be quantified by means of the radiation time-scale for the prominence 
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Fig. 7. Modulus of the eigenfunction T\ (in arbitrary units) versus the dimension- 
less distance to the slab axis corresponding to the fundamental internal fast mode 
for k z x p ~ 0.3 if a) all non-adiabatic mechanisms are considered, and b) without 
prominence radiation. The shaded region shows the location of the prominence slab. 



44i plasma flDe Moortel k Hoodl . 12004 ). 



(7-l)p^ p - 



(21) 



443 Considering fixed equilibrium parameters, the value of r r changes for different 

444 optical thicknesses of the prominence material (see regimes listed in Table I 

445 of Paper I). For Prominence (1) parameters (optically thin plasma), r r « 

446 309 s, whereas for Prominence (2) and Prominence (3) regimes (optically thick 

447 and very thick plasma), r r « 2,876 s and r r « 47,822 s, respectively, and so 

448 prominence radiation is less efficient. Obviously, r r — > oo if the radiative term 

449 is omitted. The coronal plasma is always taken optically thin. Figure [8] shows 

450 the damping time of the fundamental hybrid and internal fast modes as a 

451 function of k z x p for the different prominence optical regimes. We see that 

452 the larger the optical thickness, the larger the damping time. This effect is 

453 especially relevant at the coupling points with the external slow modes, where 

454 thermal instabilities appear if radiative losses are completely inhibited. 
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Fig. 8. Damping time versus k z x p for the fundamental a) hybrid and b) internal 
fast modes. The linestyles represent different optical thicknesses for the prominence 
plasma: Prominence (1) in solid line (this corresponds to the solid lines in Fig.[U]i and 
e), Prominence (2) in dotted line and Prominence (3) in dot-dashed line. The dashed 
line corresponds to the results when the prominence radiation is omitted (dashed 
lines in Fig. [6]i and e). The boundary conditions considered are v x = v z = T[ = 0. 

455 3. 6 Exploring the parameter space 



In this Section we investigate how the attenuation of oscillations is affected by 
changing the equilibrium parameters. The motivation of this study is based 
on the fact that the estimated values for prominence plasma parameters, 
such as temperature, density, magnetic field strength or optical thickness, 
varies from one promin ence to another, sometimes in a significant way (e.g. 
Patsourakos Viall . |2002| ). Thus, it is important for our investigation to ascer- 

462 tain the sensitivity of the damping time to the equilibrium parameters around 

463 the values considered in our previous calculations. 

464 First, we plot in Fig. [9] the ratio of the damping time to the period cor- 

465 responding to the fundamental modes as a function of equilibrium physical 

466 conditions, namely the prominence temperature, the prominence density, the 

467 magnetic field strength and the coronal temperature. The following ranges of 

468 values have been considered: 5000 K < T p < 15,000 K; 10~ n kg m -3 < p p < 

469 10- 10 kg m" 3 ; 1 G < B < 15 G; and 800,000 K < T c < 2,000,000 K. 

470 At first sight, we notice that the attenuation of fast modes is much more 
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471 sensitive to the equilibrium conditions than the damping of slow modes. The 

472 attenuation of slow modes does not change in a significant way if the equi- 

473 librium physical conditions are modified, since the obtained r D /P are always 

474 small and of the same order of magnitude. On the contrary, fast modes are 

475 highly sensitive especially to the prominence density and the magnetic field. It 

476 is noticeable that small values of r^/P are obtained for the fast modes when 

477 large densities and weak magnetic fields are considered. If the magnetic field 

478 strength is increased or the prominence density is reduced, then r D /P grows 

479 dramatically. Additionally, fast modes are again strongly affected by the cou- 

480 plings with slow modes, a fact that shows up in the form of very localised 

481 increases and decreases of r^/P. 
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Fig. 9. Ratio of the damping time to the period for the fundamental oscillatory 
modes as function of, from the left to the right, the prominence temperature, 
the prominence density, the magnetic field strength and the coronal tempera- 
ture. Computations performed considering k z x p = 1 and the boundary conditions 
v x = v z = T\ = at x = ±x c . 

482 On the other hand, we have studied the effect of considering a different heating 

483 scenario on the wave attenuation. In agreement with previous investigations 
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484 ( iCarbonell et all 12004 ; iTerradas et all 120051 ). results do not show significant 
4B5 discrepancies if different heating mechanisms are assumed. 



4B6 Finally, we have also varied the length of magnetic field lines (by modifying 

4B7 the value of x c ) and the prominence half-width, x p , in order to assess their 

488 effect on the damping time. For realistic values of both x c and x p , no significant 

489 influences appear in the results with respect to those previously discussed. It is 

490 worth to mention that prominence conduction becomes a relevant mechanism 

491 for very a small, unrealistic prominence half-width (x p < 10 km), and coronal 

492 radiation is only important for very large, and again unrealistic, length of 

493 magnetic field lines (x c > 10 6 km). 



494 3.7 Comparison with Terradas et al. (2005) 



The final check of the importance of the coronal med ium comes from the co m- 
parison between our results and those obtained by ITerradas et al.l (120051 ) in 
the case of an isolated prominence slab (see Fig. [10]). Obviously, this com- 
parison can only be performed for internal modes, since external and hybrid 
modes are not suppo rted by an isol a ted s lab. The boundary conditions as- 
sumed in the work of ITerradas et al.l (120051 ) are v x — v z — T[ — 0. According 
to the arguments given in Sect. 13.41 this condition for the perturbation to 
the temperature is the most suitable since thermal conduction is negligible 
in prominences. However, the line-tying condition at the edges of the promi- 
nence slab seems not to be the most appropriate election in the light of the 
eigenfunctions plotted in Fig. [31 Hence, our results point out that the interface 
between the prominence slab and the corona does not act as a rigid wall, and 
perturbations can be important in the corona even for internal modes. 



508 Contrary to what was shown in Paper I, in which only the fast mode was 

509 affected by the corona, in the present case both slow and fast modes of the 

510 isolated slab differ from those of a prominence-corona equilibrium. A decre- 

511 ment of the damping time is obtained for both waves in comparison with 

512 the solution of an isolated slab. The slow mode is less affected by the pres- 

513 ence of the corona but the fast mode damping time is reduced by an order of 

514 magnitude, although it is still far from the observed values. As in the longi- 

515 tudinal magnetic field case, the consideration of the corona is of paramount 

516 importance for a correct description of the behaviour of oscillations and their 

517 attenuation, although its effect on the damping of fast modes is less noticeable 
sis than in the longitudinal case. 
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Fig. 10. Period (left), damping time (centre) and ratio of the damping time to the 
period (right) versus k z for the internal fundamental slow (upper panels) and fast 
(bottom panels) oscillatory modes. The solid lines are the solutions of a promi- 
nence plus coro na equilibrium wherea s the dotted lines represent the solutions of 



an isolated slab (jTerradas et al.l . 120051 ) . 



519 4 Conclusions 



520 In this paper we have studied the wave attenuation in a system representing a 

521 quiescent solar prominence embedded in the coronal medium. The prominence 

522 has been modelled as a homogeneous plasma slab surrounded by a homoge- 

523 neous medium with coronal conditions. Magnetic field lines have been assumed 

524 transverse to the prominence slab axis and the whole system has been bounded 

525 in the field direction by two photospheric rigid walls, in order to establish a 

526 realistic length for the field lines. The attenuation of the normal modes of such 

527 equilibrium has been investigated by considering parallel thermal conduction, 

528 radiative losses and plasma heating as non-adiabatic mechanisms, and focus- 

529 ing our study on the fundamental oscillatory modes. The main conclusions of 

530 this work are summarised next. 



531 (1) Slow modes are strongly attenuated by non-adiabatic mechanisms, their 

532 damping times being of the order of the corresponding periods. Fast 

533 modes are less affected and present greater damping times. 

534 (2) The most relevant damping mechanisms are prominence radiation and 

535 coronal thermal conduction. The first one dominates the damping of in- 

536 ternal modes, while the second one is responsible for the attenuation of 

537 external modes. The combined effect of both mechanisms governs the 

538 damping of hybrid modes. Neither prominence conduction nor coronal 

539 radiation become of importance for realistic values of the length of mag- 

540 netic field lines and the prominence width. 

541 (3) The attenuation of slow modes is not affected by the value of the free 
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542 component of the wavenumber, k z . On the contrary, the behaviour of 

543 fast modes is strongly dependent on k z . 

544 (4) Thermal conduction allows energy transfer between the prominence slab 

545 and the coronal medium. Prominence radiation has an essential role in 

546 dissipating the extra heat injected from the corona and stabilises the 

547 oscillations. Thermal instabilities appear if the radiative losses from the 

548 prominence plasma are omitted or significantly reduced (e.g. caused by 

549 an increase of the optical thickness) since the plasma cannot dissipate 

550 the extra injected heat in an efficient way. 

551 (5) The damping time of fast modes is strongly sensitive to the equilibrium 

552 physical parameters while slow waves are less affected by the variation of 

553 the equilibrium conditions. 

554 (6) The presence of the corona produces a decrement of the damping time 

555 of internal modes with respect to the solutions supported by an isolated 

556 prominence slab. Nevertheless, this effect is not enough to obtain damping 

557 times of the order of the period in the case of fast modes. 



558 Considering the equilibrium parameters of Paper I, the efficiency of non- 
559 adiabatic mechanisms on the damping of fast modes is smaller in the present 

560 case. This fact suggests that the orientation of magnetic field lines with re- 

561 spect to the slab axis has a relevant influence on the attenuation of fast modes, 

562 the configuration of Paper I and the present one being limit cases. Moreover, 

563 fast modes are strongly sensitive to the equilibrium physical conditions, and 

564 it is possible to obtain small values of the damping time by considering ex- 

565 treme equilibrium parameters, such as very weak magnetic fields and very large 

566 prominence densities. In this way, fast modes show a wide range of theoretical 

567 damping times. On the other hand slow modes are always efficiently attenu- 

568 ated, with damping times of the order of their periods. This result suggests 

569 that the attenuation of prominence fast waves may be caused by other damp- 
ing mech anisms not considere d here. Some candidates c ould be resonant ab - 
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574 
575 



577 



571 sorption (lArregui et al.l . l2008l ) and ion- neutral collisions (IForteza et all 120071 ). 

572 Among these mechanisms, resonant absortion may be a very efficient damp- 

573 ing mechanism if non-uniform equilibria are considered, e.g. models with a 
transition region between the prominence and the corona. Other effects, as 
wave leakage, might only play a minor role in the damping of disturbances. 

576 Finally, future studies should take into account the prominence fine structure 
on the basis that small-amplitude oscillations are of local nature. Therefore, 



578 the investigation of the damping of fibril oscillations should be the next step. 
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